Spontaneous Spin Polarized Currents in Superconductor-Ferromagnetic Metal 

Heterostructures 
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We study a simple microscopic model for thin, ferromagnetic, metallic layers on semi-infinite bulk 
superconductor. We find that for certain values of the exchange spliting, on the ferromagnetic side, 
the ground states of such structures feature spontaneously induced spin polarized currents. Using a 
mean-field theory, which is selfconsistent with respect to the pairing amplitude x, spin polarization 
rh and the spontaneous current js, we show that not only there are Andreev bound states in the 
ferromagnet but when their energies En are near zero they support spontaneous currents parallel to 
the ferromagnetic-superconducting interface. Moreover, we demonstrate that the spin-polarization 
of these currents depends sensitively on the band filling. 
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PACS numbers: 72.25.-b, 74.50.+r, 75.75.+a 

Recently it has become possible to fabricate high qual- 
ity interfaces between superconductors {SC) and a metal- 
lic Ferromagnets (FM) . Evidently the proximity effect 
in such SC/FM hybrid structures is of scientific interest 
since it facilitates the study of coexistance between mag- 
netism and superconductivity ||]. Furthermore, it may 
also become technologically important in connection with 
magnetoelectronics and quantum computing In 
this letter we report on our study of this intriguing phe- 
nomenon in a ferromagnetic layers on a bulk supercon- 
ductor with a particular focus on one of its novel physical 
features, a spontaneously induced current in the ground 
state, which is relavent form both points of view. 

Clearly, in the present context the 'proximity effect' 
means both the leakage of superconductivity into the 
non-superconducting, ferromagnetic metal and the spin 
polarization of the superconductor near the interface. 
In the analogous case where the normal metal is non- 
magnetic (NM) this effect has been studied for a long 
time and is, by now, well understood By contrast, 
the experimental and theoretical interset in the SC /FM 
heterostructures and interfaces is more recent and the 
subject is correspondingly less well developed. Neverthe- 
less, a number of the new phenomena, associated with 
Cooper pairs in an exchange field, have been indetified. 
For instance, it has been found that, as opposed to the 
SC/NM case, the pairing amplitude x does not de- 
cay exponentially to zero on the ferromagnetic side of 
a SC/FM interface but oscilates, with a slowly decreas- 
ing amplitude, as a function of the distance d from the 
interface (|]-||]. These oscillations turn out to be mani- 
festations of the effect first studied by Fulde and Ferrel 
pO| and Larkin and Ovchinikov [|ll], often refered to as 
FFLO, and some of their more striking consequences, 
such as the oscillation of the superconducting transition 
temperature Tc as a function of the thickness of the FM 
layer has been observed experimentally p2[. 



Our work is particularly relavent to the predicted 
|l3| , [T^ and observed jTs) Andereev bound states in clean 
thin ferromagnetic films sandwiched between two super- 
conductors. Interestingly, these states were found to form 
part of the ground or equilibrium states and to give rise 
to so called 7r-states of the SC / FM / SC junction. One 
of the aim of this letter is to argue that a hitherto over- 
looked salient feature of such FFLO 7r-jucntion is a spon- 
taneous current parallel to the FM/SC interface. We 
shall also show that the spin polarization of such current 
depends sensitively on whether the state has, or has not, 
particle-hole symmetry. 

To investigate the occurence spontaneous currents and 
their polarization in a SC /FM heterostructure we need 
a model simple enough to be solved, at least in a mean- 
field approximation, for the magnetization m, pairing 
amplitude x and spontenous current js self-consistently. 
Hence, for the purpose at hand, we have adopted a sin- 
gle orbital, nearest neighbour hopping, negative U Hub- 
bard model for describing the semi-infinite superconduc- 
tor and continued the same Hamiltonian into the ferro- 
magnetic layer with the U set equal to zero and the site 
energies Sic exchange split. Moreover, we consider the 
simplest geometry, depicted in Fig.^ where a magnetic 
field in one direction, the vector potential and a current 
in another and a spatial modulation in a third orthogo- 
nal direction can be realised. Hopefully, while simplify- 
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FIG. 1: Schematic view of the superconductor-ferromagnet 
interface. Directions of the magnetic field (_B) as well as vector 
potential (^4) and current (J) are indicated. 
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ing the calculation this effectively 2D system will have 
much in common with its iD counterpart such as a layer 
of ferromagnetic metal deposited on a superconducting 
substrate |l^ and corresponding sandwich structures 

In short, our model Hamiltonian is given by 

H = X^t^^J' + (^^'^ ^ M)^y]c+ Cj> + Y ^fii^ni-y (1) 

ijcr ia 

where, in the presence of a vector potential A{r), tij = 

—te ^'^ '^^'^^ '^^ for nearst neighbour lattice sites at spa- 
tial positions Ri and Rj , the site energies Eia are on the 
superconducting side and equal to ^E^^a on the ferro- 
magnetic side, is the chemical potential, Ui is Us < 
in the superconductor and zero elswhere, c^, {cia) are 
the usual electron creation (annihilation) operators and 
flier — c^^Cicr- Notc that the above description of the elec- 
trons with charges e includes a coupling to a magnetic 
field B{r) = V x A(r). Evidently such provision will be 
necessary for calculating the effects of currents on the 
electronic states. 

In what follows we shall study the above model in 
the Spin-Polarized-Hartree-Fock-Gorkov (SPHFG) ap- 
proximation. We shall work in the Landau gauge where 
B = (0,0, J5^(a:)) and hence A= (0, ^j;(a;), 0). Further- 
more, we assume that the effective SPHFG Hamiltonian 
is periodic in the direction parallel to the inteface and 
therefore we work in k space in the y direction but in 
real space in the x-direction (see Fig.0). Labeling the 
planes in Fig.|l| by integers n and m at each ky point of 
the Brillouin zone we shall solve the following SPHFG 
Nambu, spin (a) and layer index (n) matrix equation: 

E H:Z,{u;,ky)Gl',^{cj,ky) = 6n^S^f, (2) 

rn' ,^,ky 

where the only non-zero elements are: H^ln ^^'^ ^nm = 
{u - \iyEex ± /i ± tcos{ky =F eA{n)))Snm ± tSn,n+i for 
the upper and lower sign respectively, = H^l^ and 

H^'^ — with a replaced by — cr in both cases, 

fjl2 P721 _ f734 _ ff43 A A anH /^'^P 

is corresponding retarded Green's function {GF). 
As usual, selfconsistency is assured by the relation: 

ky 

-[/„E / du;hmGll{LO,ky)fiLo) (3) 

ky 

where /(w) is the Fermi distribution function. Moreover, 
the FM order parameter is given by 

mn = ^{nn^ - n„x) ^ x 

ky 

J cLolmiGliico, ky) - G33„(a;, ky))f{io) (4) 



Since our model includes a vector potential the solution 
of Eq.^ will imply a current. For spin up electrons, in the 
y-direction this can be calculated from the relation: 

Jj,l(j)(n) = — 2et E] sin{ky — eAy{n)) x 

ky 

J dujhmGlli''\u;,ky)f{co) (5) 

which follows from the continuity equation for the charge. 

Finally, the above current will give rise to a vector 
potential Anewi^) which will have to be used to update 
A{'r) in Eq.^ at the end of each selfconsistency cycle. We 
calulated this new vector potential by solving numerically 

Ampere's law, j^i'^"' = ^ 47rJj^(x), which for the lattice 
problem at hand, is 

Ay{n + 1) - 2Ay{n) + Ay{n - 1) = -'iTrJy{Ti) (6) 

We have solved Eqs using an appropriately simpli- 
fied 'principal layers method' ||l7| which we shall describe 
elswhere ||T^. The rest of this letter is a brief summary 
of our resuls. 

Firstly, since we determined the order parameters, 
Xn and m„, on both sides of the inerface fully self- 
consistently, we were able to study both a superconduct- 
ing and a magnetic proximity effects. Although A„ — 
on the ferromagnetic side, due to the fact that U — 0, 
the pairing amplitude x„ = (c„|C„f) does not have to 
be, and indeed it turns out not to be, zero. Given the 
well understood effect of the exchange field in a bulk 
superconductor [|l0[ |ll| it is not all together surprising 
that we find that on entering into the ferromagnet \n 
oscillates as a function of the distance from the inter- 
face. In fact our numerical results fit the analytic formula 
xix) ~ sin{x / ^f) / {x / £_f) , where = t/E^x is the fer- 
romagnetic coherence length, and hence are fully consisi- 
tent with those of Ref . ^, ^ . As it is well known 

at a NM /SC inerface Xn decays exponentially in the nor- 
mal region with the decay length ^5. Thus, although the 
supercoducting coherence length ^5 >> due to the 
functional form, polynomial as opposed to exponential, 
of the decay Xn enters more deeply into the ferromag- 
net then into the normal metal. An other novelty of the 
FM / SC interface, when compared with its well studied 
NM/SG counterpart, is the entry of magnetism into the 
superconductor. We find that the spin polarization does 
not oscillate but decays monotonically on a length scale 
of order of the ^p- As shown in Figj^, somewhat supris- 
ingly, m„ oscillates on the ferromagnetic side but dacays 
exponentially on the superconducting side of our hybrid 
structure. This behavior is in full agreement with the 
results of Demler et al. |^ . 

Of course, the most remarkable feature of the above 
solution is that the iterations of the SPHFG equations 
frequently converge to a finite value of the current jy (n) 
even though the external vector potential is zero. We 
have checked that the exsistence of such spontaneous cur- 
rent lowers the energy of the system. As shown in FigJ3, 
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FIG. 2: Magnetization and the pairing amplitude as a func- 
tion of the distance from the interface for Eex = 1.273 and 
Us ~ —2, which gives As = 0.376 in units of the hopping 
integral t. Note that we plot — Xn which corresponds to posi- 
tive An for Un > 0. Inset: Example of the semiclassical paths 
corresponding to Andreev bound states. 



it flows in the positive y direction on the ferromagnetic 
side, and in negative in the superconductor and sums, 
reassuringly, to zero over all layers. We have also found 
numerically, that there is a magnetic flux $ w *i>o/2, <&o> 
being the flux quantum, associated with this current dis- 
tribution. 
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FIG. 3: Total current Jy°'(n) = JyT{n) + Jyi{n) and the total 
density of states at the Fermi level Pnisp)- For comparison 
the DOS for Eex = {pZ{£f)) is also shown. Inset: Inte- 
grated over the FM side total DOS pfm{ijj)- Dashed (solid) 
line corresponds to the solution with (without) spontaneous 
current. The parameters are the same as in the Fig.^ 

To shed light on the origin of the above spontaneous 
current, in Fig.^ we also show the layer resolved den- 
sity of states, namely the spectral function Pniep) = 
Efc„ Im(Gii„(eF, ky) + Gf^{eF, ky)), at the Fermi en- 
ergy ep- Clearly, the oscillations of the layer resolved cur- 
rent tracks those of the density of states. Further insight 
follows if we compare this with the corresponing quantity, 
denoted by p° (^i;-) and represented by the dot line, in the 
case where the exchange spliting, E^x, is set equal to zero. 
Since the rise and fall of p'^{ep) across the ferromagnetic 
layer can be readily inetrpreted as the order parameter 



amplitude of an Andreev bound state, corresponding to 
the semiclassical path depicted in the inset of the Fig||, 
we can regard Xn in Fig.^ and pnisp) in Fig.^ as an indi- 
cation that a similar bound state is formed in the much 
more complicated case of finite exchange field. In fact 
we can indentify such i^FLO-Andreev bound states as 
peaks in the full quasiparticle density of states for the 
ferromagnetic layers. As might be expected these are 
exchange split and move around as a function of the ex- 
change field Ef-x- Investigating the correlation between 
such bound states and the current carrying capacity of 
a solution we find that in the ground state a current 
flows only when the energy of one of the FFLO- Andreev 
bound state is near the Fermi energy. Furthermore, in 
the presence of the current, this state splits, thus lower- 
ing the total energy of the system. An example of such 
zero energy bound state is depicted in the inset of the 
Fig.|. 

Zero energy Andeev bound states in SC / N M / SC juc- 
tions usually lead to a so called 7r-states of the two super- 
conductor in which the phases of their order parameters 
differ by tt. Recently, Chtchelkatchev et al. |^ sug- 
gested that the same is true for a SC / FM / SC junction 
in the presence of fully developed FFLO phenomena. In- 
terestingly, although the structure we have been studying 
has only one superconducting region it turns out to dis- 
play properties analougous to those of such vr-junctions. 
To see these we studied the order parameter x_g at the 
surface of the ferromagnetic layer opposite to the super- 
conductor. Since Xn is negative in superconductor when 
X-9 is positive we may describe the system as being in a 
TT-state. In FigJJ we display our results for x-9 as a func- 
tion of the dimensionless exchange held 8 = SdEex/T^t 
{d is the number of FM layers). Evidently, there are 
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FIG. 4: Pairing amplitude at the surface of the FM layer 
oposite to superconductor (n — —9) as a function of the di- 
mensionless parameter O = SdE^x/TTt. Note pinning of the x 
to zero as long as current flows. 

regions of for which the sytem can be said to be in 
TT-state, in close agreement with the analougous results 
of Chtchelkatchev et al. 

Remarkably, the states with spontaneous currents cor- 
respond to regions of the exchange held Q where the or- 
der parameter x-9 is near zero. As can be easily read 
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off from these plots, in these regions the presence or ab- 
sence of the currents have a dramatic effect on the order 
parameter. Evidently, the spontaneous current pins the 
order parameter on the non-superconducting side of the 
ferromagnetic layer (x-g) to zero. This appears to stab- 
lise the zero energy i^i^iO- Andreev bound state. Given 
their physical origin one might expect the above sponte- 
neous currents to be spin polarised. This is indeed the 
case. In fact, we find that the degree of spin polariza- 
tion is largely detemined by the difference in the spin 
up and spin down densities of state at the Fermi energy: 
Apn = pn^isp) — Pniisp)- In the above calculations we 
have assumed a half filled band, that is to say particle 
hole symetry, and hence we have found no difference be- 
tween the spin up and spin down currents. However, 
further calculations, away from particle hole symetry, re- 
vealed much larger spin polarizations. An example of 
this is reported in Fig||. 




FIG. 5: Spin polarized current AJy{n) = Jyi{n) — Jyi{n) 
and the difference of the spin up and spin down density of 
states at the Fermi level Apn{£F) ~ p„|(eF) — pniisp) as a 
function of the distance from the interface. Here the band 
filling is 0.636. 



In summary, we have demonstrated that zero energy 
Andreev bound states in ferromagnetic layers deposited 
on a superconducting substrate will carry currents in 
the ground state of such a hybrid structure. Moreover, 
we found that such states will form only in certain re- 
gions of the exchange field E^x and layer thicknes d 
phase diagram. In particular, we investigated the cases 
< Eex/t < 3, < d < 20. Finally we found that the 
spin polarization of the current is closely related to the 
difference in the spin up and spin down density of states 
at the Fermi level. 
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